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We report an experimental study of NaCu202, a Mott insulator containing chains of edge-sharing 
Cu04 plaquettes, by polarized x-ray absorption spectroscopy (XAS), resonant magnetic x-ray scat- 
tering (RMXS), magnetic susceptibility, and pyroelectric current measurements. The XAS data 
show that the valence holes reside exclusively on the Cu^^ sites within the copper-oxide spin chains 
and populate a d-orbital polarized within the Cu04 plaquettes. The RMXS measurements confirm 
the presence of incommensurate magnetic order below a Neel temperature of Tjv — 11.5 K, which 
was previously inferred from neutron powder diffraction and nuclear magnetic resonance data. In 
conjunction with the magnetic susceptibility and XAS data, they also demonstrate a new "orbital" 
selection rule for RMXS that is of general relevance for magnetic structure determinations by this 
technique. Dielectric property measurements reveal the absence of significant ferroelectric polar- 
ization below Tjv, which is in striking contrast to corresponding observations on the isostructural 
compound LiCu2 02. The results are discussed in the context of current theories of multiferroicity. 



PACS numbers: 75.25.-j, 75.85.+t, 71.27.+a, 75.50.-y 

I. INTRODUCTION 

Stimulated in part by the discovery of high- 
temperature superconductivity, copper oxides with spin- 
1/2 networks ranging from zero to three dimensions have 
been the focus of a large amount of research during 
the past two decades. Among these, Mott-insulating 
compounds that comprise spin chains with edge-sharing 
Cu04 plaquettes have been of particular recent interest 
as models of frustrated magnetism in one dimension. In 
these materials, the Cu-O-Cu bond angle is close to 90°, 
so that the magnitude of the superexchange coupling be- 
tween nearest-neighbor Cu^^ spin-1/2 ions, Ji, is smaller 
than that of the next-nearest neighbor coupling, J2.— De- 
pending on the sign and magnitude of the ratio J1/J2, as 
well as on the relative magnitudes of the spin-lattice cou- 
pling and the interchain exchange interactions, a variety 
of ground states have been observed in different mem- 
bers of this family, including spiral magnetism^i^ii^i^ 
collinear antiferromagnetismj^ and spin-Peierls order 
The quasi-one-dimensional spin correlations have been 
observed to profoundly affect charge excitations across 
the Mott-Hubbard gapi^S In addition, chargeiiJ^ and a 
spin-density modulation-^- have been reported in materi- 
als with doped edge-sharing copper-oxide chains. 

The recent discovery of ferroelectricity associated with 
helical magnetic order ( "multiferroicity" ) in the insulat- 



ing edge-sharing copper-oxide chain compounds LiCu202 
(Ref. M) and LiCuV04 (Ref. HI) has stimulated another 
wave of research on this family of materials, and it has 
sparked a controversy about the origin of this effect. The 
magneto-electric coupling is generally consistent with 
models according to which multiferroicity is an intrinsic 
consequence of spiral magnetism i^^i^^i^^i^^ While some 
aspects of the experimental data on both compounds 
agree with specific predictions of these models, ^'''^''•^^ dis- 
agreements with other predictions have also been pointed 
outi ^^^i^^i^^ An alternative explanation invokes orbital po- 
larization of defects generated by inter-substitution of Li 
and Cu (whose ionic radii are closely similar) as the ori- 
gin of the electric polarizatioui^^iSiiS This scenario is 
supported by the observation of substantial deviations 
from stoichiometry and substitutional disorder in both 
LiCu202 and LiCuV044'iS2, However, the applicability of 
this model to LiCu202 has been disputed based on the 
results of a recent study on single crystals with nearly 
stoichiometric compositioni^i^ 

We report the results of an experimental investiga- 
tion of single crystals of NaCu202, a Mott insulator 
that is isostructural and isoelectronic to LiCu2 02 (Fig. 
1). Like LiCu202 and LiCuV04,-''''^ NaCu202 exhibits 
helical magnetic long-range order at low temperatures 
due to competing ferromagnetic nearest-neighbor and 
antiferromagnetic next-nearest neighbor superexchange 



2 




FIG. 1: (Color online) Sketch of the orthorhombic crystal 
structure of NaCu2 02 (space group Pnvna). The room- 
temperature lattice parameters are a—&. 2087 A, 6=2.9343 A 
and c=13.0548 A.^^ 

interactions.—'^ The main difference to the Li-based com- 
pounds is the large difference between the ionic radii of 
Na and Cu, which disfavors inter-substitution of both 
elements. •i'^ Experiments on NaCu2 02 thus offer the 
chance to elucidate the intrinsic properties of this class 
of compounds, without interference from disorder. 

This paper is organized as follows. Section II de- 
scribes polarization-dependent x-ray absorption spec- 
troscopy (XAS) experiments designed to clarify the or- 
bital character of the valence electrons in NaCu202. Sec- 
tion III contains the results of susceptibility and resonant 
magnetic x-ray scattering (RMXS) experiments aimed at 
elucidating the magnetic properties of NaCu202. Section 
IV summarizes measurements of the electrical polariza- 
tion and the dielectric constant. In the concluding Sec- 
tion V, the results are contrasted with observations on 
LiCu202 and discussed in the light of current theories of 
multiferroicity. 



II. ELECTRONIC STRUCTURE 

We begin by describing XAS measurements near the 
Cu L-absorption edge aimed at determining the orbital 
character of the valence-electron states in NaCu202. The 
NaCu2 02 lattice structure contains copper ions in two 
non-equivalent sites [Cu(l) and Cu(2) in Fig. 1]. The 
chemical coordination and bond lengths of these ions in- 
dicate valence states of and 2+ (with electron con- 
figurations 3(i^° and 3d^, and spins and 1/2), respec- 
tively. However, a recent XAS study of the isostructural 
LiCu202 had come to the conclusion that the Cu(l) ions 



contain an intrinsic density of holes, independent of Li- 
Cu inter-substitution.— We have checked the validity of 
this conclusion on NaCu2 02, a compound in which the 
latter complication is manifestly absent. 

The experiments were performed on an untwinned, 
plate-like single-crystal with surface area 3x2 mm'^ and 
thickness 100 /im, which had been prepared by a self-flux 
technique already described elsewherci^ The sample sur- 
face is parallel to the ah plane of the orthorhombic crystal 
structure and hence includes the CUO4 plaquettes (Fig. 
1). The XAS data were taken at beamline 4-ID-C of 
the Advanced Photon Source in total-fluorescence-yield 
mode at a temperature of 150 K. The beam polarization 
was horizontal, and the sample was rotated such that 
data could be taken with different angles, 0, between the 
incident wave vector ki and the sample surface (Fig. [2]). 
All spectra were normalized to the same intensity 10 eV 
above the Cu L2 absorption edge. 

Figure 2 shows a series of XAS spectra obtained for 
different 6. The features in the energy ranges 930-938 eV 
and 950-958 eV correspond to the CU-L3 and L2 absorp- 
tion edges, respectively. At the L3 edge, three peaks are 
observed. The energy splitting of 2.7 eV between the 
first peak (Pi) at 930.5 eV and the second peak (P2) at 
933.2 eV allows us to attribute Pi to the Cu^+ ions on 
site Cu(2), and P2 to the Cui+ ions on site Cu(l),^ The 
additional peak (P3) at 936.7 eV is related to an excita- 
tion into non-local states. '^'^ When ki is turned from in- 
plane to out-of-plane, the intensity of Pi is progressively 
enhanced, indicating that the valence hole on Cu^+ oc- 
cupies c?-orbitals polarized within the CUO4 plaquettes. 
The same polarization dependence is observed at the Cu 
L2 absorption edge. These findings agree with density 
functional calculations for the structurally nearly iden- 
tical LiCu202 compound, which predict that the states 
near the Fermi level are hybrids of the Cu d^2_y2 orbital 
and O Px.y orbitals in the Cu02 chainsJ^i^^ 

The intensity of peak P2 originating from the Cu(l) 
ions is independent of 9 within the experimental error 
(inset in Fig. 2), as expected for the spherically symmet- 
ric 2p^3d^° ground-state and 2p^3d^°4s^ final-state con- 
figurations in the XAS excitation of Cu^"*" . We note that 
our data are quite different from corresponding data on 
LiCu202 that revealed a strong polarization dependence 
of the P2 peak intensityf22, Based in part on electronic- 
structure calculations, the authors of Ref. [l^ had at- 
tributed this observation to an intrinsic density of holes 
on the nominally monovalent copper ions, independent of 
the deviations from stoichiometry in LiCu202. Our data 
on the stoichiometric, but structurally nearly identical 
NaCu202 do not support this interpretation. Rather, 
they indicate that the Cu d-orbitals on the Cu(l) sites of 
both compounds are fully occupied, as expected for the 
1+ valence state, and that the polarization dependence 
of the corresponding XAS peak in LiCu202 originates 
from defects created by Li-Cu inter-substitution. 



3 



2.5 




O '' ' ■ ■ ■ - ' 

920 930 940 950 960 

Energy (eV) 



FIG. 2: (Color online) Polarization-dependent x-ray absorp- 
tion spectra around the Cu L3 and L2 absorption edges of 
NaCu202. 9 denotes the angle between the incident wave vec- 
tor and the sample surface, which was parallel to the ab plane 
(see lower right inset). The origins of the peaks Pi, P2, P3 
are discussed in the text. The inset shows the 6'-dependence 
of the intensity of peak P2 on an larger scale. 



III. MAGNETIC PROPERTIES 

Figure 3 shows the magnetic susceptibility, Xj mea- 
sured on an untwinned crystal from the same batch as 
the one described above. The broad maximum around 
T=50 K is characteristic of short-range antiferromagnetic 
spin correlations in low-dimensional magnets. Upon cool- 
ing below the Neel temperature T;v=11.5 K, x exhibits 
a sharp drop for magnetic fields, H, along h or c, while a 
slight upturn is observed for H || a. This implies that the 
direction of the magnetic moment in the Neel state is pre- 
dominantly in the he plane, consistent with a 6c-polarized 
spiral inferred from neutron powder diffraction and nu- 
clear magnetic resonance (NMR) measurements A 
previously reportect^i^. anomaly in x{T) for T 8 K was 
not observed in our samples. Measurements of x(H) up 
to 7 T (not shown) did not reveal any anomalies indica- 
tive of spin- flop transitions. 

RMXS measurements with photon energies near the 
Cu L3 edge were performed in order to further investigate 
the spin correlations below Tjy. The experiments were 
performed at the U46-PGM1 beamline at the BESSY 
II synchrotron in Berlin, Germany, using a two-circle 
diffractometer with horizontal scattering geometry de- 
signed at the Freie Universitat Berlin. Silver paste was 
used to glue the crystal onto a copper goniometer at- 
tached to a continuous flow cryostat, such that the a- and 
&-axes were parallel to the diffraction plane. A 800 nm 
thick aluminum foil provided sufficient heat shielding to 
reach sample temperatures as low as 3 K. The incident 
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FIG. 3: (Color online) Temperature dependence of the mag- 
netic susceptibility along the principal crystallographic axes 
measured with a magnetic field H=1000 Oe. The inset high- 
lights the data at low temperature. 

beam was polarized either horizontally (H) or vertically 
(V), i.e. parallel or perpendicular to the scattering plane. 
The diffracted beam was detected without polarization 
analysis. 

At temperatures below T/y, a resonant diffraction peak 
was observed at wave vector transfer Q = (0.5 ^ 0), with 
incommensurate component ^ 0.227 along the spin 
chains. (The components of Q are indexed based on the 
orthorhombic crystal structure; see Fig. 1.) The origin of 
the peak was identified as magnetic based on the temper- 
ature dependence of its intensity (Fig. 4a). The RMXS 
data thus confirm the propagation vector extracted from 
powder neutron diffraction measurements The width of 
the diffraction peak, measured along the propagation vec- 
tor, is temperature independent and yields a lower bound 
of 3500 A on the magnetic domain size in the ab plane. 
The temperature variation of the incommensurate wave 
vector below T/v (Fig. 4b) is much smaller than the corre- 
sponding variation in LiCu2 02, where a two-step transi- 
tion to magnetic long-range order was observed,^i25 The 
energy dependence of the (0.5 ^ 0) peak intensity (Fig. 
5) shows two minima at 931.6 and 934.3 eV, which can 
be attributed to strong absorption at the L3 edges of 
Cu^+ and Cu"*" (dotted lines in Fig. 5). Similar data 
were reported for LiCu202 (Ref. [s^. The resonant mag- 
netic Bragg reflection was also observed when the pho- 
ton energy was tuned to the Cu L2-absorption edge of 
NaCu202, but an equivalent experiment at the Cu K- 
edge yielded a null result. This is consistent with the ex- 
pectation that the resonant enhancement of the RMXS 
cross section at transition-metal L-edges is much larger 
than at K -edges. 

We now discuss the dependence of the RMXS cross 
section on the incident beam polarization ci and on the 
aziniuthal angle The scattering geometry and the def- 
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FIG. 4: (Color online) (a) Temperature dependence of the 
intensity of the magnetic Bragg reflection at Q=(0.5 ^ 0), 
extracted from fits to RMXS data, (b) 6-axis component of 
Q as a function of temperature. 



inition of are illustrated in the inset of Fig. 6. Note 
that Q of the magnetic Bragg reflection approximately 
bisects the angle between the a- and 6-axes, and that 
for = 0° the final wave vector is nearly along b 
(and the outgoing polarization e/ is hence nearly _L b). 
For 5* — 180°, on the other hand, k/ is nearly along 
a (if nearly _L a). According to the standard treatment 
of RMXSf22 the magnetic scattering amplitude is propor- 
tional to [(e/ X ii) -m], where m is the magnetic moment. 
The susceptibility indicates that m is confined to the be 
plane below Tn (Fig. 3). For = 0° and V (H) incident 
polarization, one would hence expect the scattering am- 
plitude to be mainly proportional to mt (rric), while for 
^ = 180° it should nearly vanish (be sensitive to both 
nib and rric). Since the susceptibility further indicates 
comparable magnitudes of nib and nic, the ratio Iv/Ih 
is expected to be of the order of unity for ^E* = 0° , and to 
nearly vanish for ^ = 180°. 

In striking contrast to this expectation, we have found 
that Iv is about three orders of magnitude smaller than 
Ih for ^ = 0° (Fig. 6a). Moreover, the azimuthal- 
angle dependence of the ratio Iv/Ih (which should not 
be strongly affected by distortions arising from sample- 
size or surface effects) is well described by the function 
Asin^vE' with A ~ 1 (Fig. 6b), which would imply 
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FIG. 5: (Color online) Energy dependence of the scattered 
intensity at the (0.5 0.227 0) magnetic reflection near the Cu 
1/3 absorption edge. The data were taken at T = 3 K. The 
dotted lines at 931.6 eV and 934.3 eV represent the energies 
of the Cu'^"'' and Cu^"*" absorption edges, which were extracted 
from the widths of longitudinal momentum scans through the 
magnetic reflection. 



vria ~ mb according to the standard description of 
RMXSi^ A recent generalization of the theory of RMXS 
that takes the symmetry of the valence electron wave 
function into account^^ resolves the apparent contradic- 
tion between the susceptibility and RMXS data. In the 
case at hand (which is explicitly treated in Ref. 38), it 
was shown that a selection rule precludes excitations of 2p 
core electrons into the partially occupied, planar Sd^^-y^ 
orbital of Cu^+ (see Section II) by vertically polarized 
photons. This implies that ly = for both = 0° and 
180°. Our data thus confirm the revised analysis of Ref. 
[38l and highlight its importance for magnetic structure 
determination by RMXS. 



IV. DIELECTRIC PROPERTIES 

Measurements of the electrical polarization P and the 
dielectric constant e were performed on a NaCu202 crys- 
tal with dimensions ?» 5 x 2 x 0.4 mm^. Due to the small 
thickness of the sample, the experiments were restricted 
to the c-axis components of both quantities. Electrodes 
were fabricated by covering the ah faces of the crystal by 
silver paste. 

The sample was loaded into a 14 T cryomagnet, and 
cooled down to 2 K while applying a poling electric field 
of 250 kV/m. Then, Pc was determined through a time- 
integral of the pyroelectric current during a temperature 
sweep rate of 20 K/min. Sc was measured with an LCR 
meter operating at a frequency of 100 kHz with a tem- 
perature sweep rate of +2 K/min. 

As shown in Fig. 7a, neither a spontaneous nor 
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FIG. 6: (Color online) (a) Longitudinal reciprocal-space scans 
through the Q=(0.5 ^ 0) reflection taken with horizontal (H) 
and vertical (V) polarization of the incident x-ray beam. 
The data for V-polarization were multiplied by 100. (b) 
Azimuthal-angle ('I') dependence of the ratio of Bragg peak 
intensities measured with V- and H-polarized incident beams. 
The line is the result of a fit to the function As'm^ '5, as dis- 
cussed in the text. A small ofi^set in *1/ was also included in 
the fit in order to account for imperfect sample alignment. 
The inset illustrates the definition of ^' and the scattering 
geometry of ^ = 0°. 



a magnetic-field-induced electric polarization were de- 
tected for magnetic fields H < 14 T within the exper- 
imental sensitivity of 0.3/iC/m^. This upper bound is 
more than an order of magnitude lower than the Pc ob- 
served in LiCu2 02 for T < The dielectric con- 
stant exhibits a subtle downturn upon cooling below 
(Fig. 7b). Effects of this magnitude have also been ob- 
served in other antiferromagnetic insulators^ and may 
be a consequence of conventional spin-phonon coupling. 
Sharp spikes in the temperature dependence of Sc such 
as the one observed at the Neel temperature of LiCu2 02 
(Refs. [l3f20l ) are not detected within the experimental 
sensitivity. 



FIG. 7: (Color online) (a) Electrical polarization and (b) di- 
electric constant measured along the c-axis for H = and 
14 T applied along the three principal crystallographic direc- 
tions. The data in panel (b) were arbitrarily shifted in vertical 
by < 5% for clarity and for taking into account small changes 
in the electrodes connection. 



V. CONCLUSIONS 

At first sight, the absence of multiferroicity in 
NaCu202 appears to contradict models that predict 
a ferroelectric polarization of the form P cx x 
(Si X Sj) (where hij is a vector connecting nearest- 
neighbor spins Sij) as an intrinsic consequence of spi- 
ral magnetismJ^iiiii^ii^ Neutron powder diffraction and 
NMR measurements have indicated a 6c-polarized spiral 
with propagation vector along which is also con- 
sistent with the susceptibility and RMXS data presented 
here. The models thus predict a ferroelectric polarization 
along c, which was not observed (Fig. 7). Refinement 
of a comprehensive set of single-crystal magnetic neu- 
tron diffraction data has, however, revealed a more com- 
plex magnetic structure in which the spin chains within 
the unit cell (Fig. 1) exhibit spiral states with oppo- 
site helicities^S so that their ferroelectric polarizations 
(if present) are expected to cancel out. Our data are 
consistent with such an antiferroelectric state below T/v, 
and hence do not contradict the theoretical work of Refs. 
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We have also shown that the Cu(l) sites of the 
NaCu202 structure exhibit a full-shell configuration and 
are hence electronically inactive. This agrees with prior 
considerations based on the chemical coordination of 
these sites, but disagrees with the conclusions of a recent 
XAS study of LiCu202 according to which they exhibit 
an intrinsic density of holesi^ The orbital polarization 
detected in the Cu(l) XAS spectra of LiCu202 is thus 
presumably a consequence of defects generated by Li-Cu 
inter-substitution. In view of the nearly identical lattice 
structures of LiCu202 and NaCu202, it seems natural to 
associate these defects with the ferroelectric polarization 
detected in LiCu202, but not in NaCu202. We cannot 
rule out, however, that the qualitative difference between 
the dielectric properties of both compounds originates in 
a subtle difference between their magnetic structures. 

Finally, the simultaneous determination of the orbital 
occupation of the valence electrons of NaCu2 02 by XAS 
and their magnetic structure by RXMS allowed us to 



experimentally confirm an "orbital" selection rule for 
RMXS recently proposed on the basis of a theoretical 
analysis. "^^ Especially for non-coUinear spin structures, it 
is thus important to consider the orbital occupation be- 
fore drawing conclusions about the spin polarization from 
RXMS data. 
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